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We present our theoretical results on the electronic state in vanadium spinel oxide AIV2O4. 
The material is a mixed- valent system with the average valence V 2 5+ , and V cations consti- 
tute a frustrated pyrochlore structure. It shows a structural transition at ~ 700 K, leading 
to the formation of seven V-sites clusters — heptamers. We study the electronic state of 
the heptamer by explicitly taking account of orbital degree of freedom as well as electron 
correlations. We show that the ground state of the heptamer for realistic parameters becomes 
spin-singlet because of strong tr-type bonding states of ti g orbitals. The temperature depen- 
dence of the magnetic susceptibility in experiments is naturally understood by this singlet 
formation in heptamers. Our results indicate that in AIV2O4 orbital physics is relevant to 
stabilize a cluster-type singlet state instead of a previously proposed charge-ordered state 
with valence skipping. 

KEYWORDS: AIV204, spinel, pyrochlore lattice, heptamer, spin singlet, exact diagonalization, 
multiorbital Hubbard model 



1. Introduction 

Geometrical frustration in strongly correlated systems has attracted much interest for 
decades. Many well-studied examples are found in frustrated antiferromagnets, where frustra- 
tion results in nearly degenerate ground-state manifolds of a large number of different spin 
configurations. 1,2 In some cases, the degeneracy remains down to the lowest temperature and 
leads to exotic phenomena such as a liquid state and a glassy state. In general, however, 
nature does not favor the degeneracy, and tries to find a way of lifting it and to select a 
unique, nondegenerate ground state. It is very intriguing to clarify the mechanism of lifting 
the degeneracy. 

The pyrochlore lattice is a typical frustrated structure in three dimensions, and the mag- 
nets on this lattice structure have attracted much attention because of their severe frustration. 
The pyrochlore lattice consists of a network of corner-sharing tetrahedra, and is regarded as a 
three-dimensional analogue of the Kagome lattice. In fact, the pyrochlore lattice can be viewed 
as an alternative stacking of Kagome and triangular planes as shown in Fig. 1. This lattice 
structure is found in many real compounds, for instance, in spinels, pyrochlore compounds 
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Fig. 1. Schematic picture of the lattice structure of AIV2O4. For simplicity, only V sites which consti- 
tute the pyrochlore structure are shown. There are three crystarographycally inequivalent V sites 
(VI, V2 and V3), and V-V bond lengths are substantially different below the structural transition 
temperature T c ~ 700 K. Trimers are shown by bold lines. Two trimers sandwiches a V2 site 
(connected by gray lines), forming a heptamer. See the text for details. 



and cubic Laves-phase intermetallic compounds. 

One of the most fundamental problems in the frustrated pyrochlore systems was studied 
by Anderson. 3 He examined the ground state and its degeneracy in the Ising model on the 
pyrochlore lattice, which mimics a possible charge-ordering phenomenon in magnetite Fe304, 
i.e., the so-called Verwey transition. 4 He pointed out that the ground state of this model does 
not exhibit any long-range order and suffers from macroscopic degeneracy. This shows that the 
frustration may give rise to nontrivial phenomena even in charge degree of freedom. Recently, 
this 'charge frustration' has also been an issue in several other mixed-valent materials such as 
CuIr 2 S 4 5 and LiV 2 4 . 6 ' 7 

A spinel oxide AIV2O4 is a typical system where this charge- frustration phenomenon is 
addressed. In this compound, vanadium cations form the pyrochlore lattice and their average 
valence is 2.5+. The vanadium cations locate at the octahedral positions, and therefore, five- 
fold d levels of vanadium are split by the crystal field into low-energy three-fold t2 g levels and 
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high-energy two-fold e g levels. Hence, we have 2.5 d electrons in the three-fold tig levels on 
average. Because the tig levels are partially filled, orbital degree of freedom is also relevant 
in this compound in addition to charge and spin degrees of freedom. Therefore, AIV2O4 is an 
interesting system to study the effects of geometrical frustration and the roles of charge, spin 
and orbital degrees of freedom. 

AIV2O4 shows a phase transition at T c ~ 700 K. At this temperature, a structural change 
is observed from the high-temperature cubic phase to the low-temperature rhombohedral 
phase accompanied by a doubling of the unit cell along the [111] direction. 8,9 This transition 
was interpreted as a charge ordering with valence skipping in which the charge distribution 
is differentiated between two different [111] planes, namely, Kagome and triangular planes 
in Fig. 1; it is supposed that V cations in Kagome layers take the valence 2+ and those in 
triangular layers take 4+. However, there remain two unresolved issues in this charge-ordering 
scenario. One is the discrepancy between the periodicity of the lattice structure and the charge 
order. The lattice unit cell is doubled in the [111] direction, which is twice as long as that 
of charge order. The other unresolved issue is the temperature dependence of the magnetic 
susceptibility. The susceptibility shows a sudden drop at T c followed by a Curie-like component 
in the low-temperature phase. It is difficult to explain this temperature dependence on the 
basis of the charge-ordering scenario. 

Recently, the lattice structure of this compound was reinvestigated. 10 It was found that 
below T c , there is an additional superlattice structure within the Kagome planes, i.e., a 'trimer' 
formation as shown in Fig. 1. The trimers are paired between neighboring Kagome planes, 
which leads to the doubling of the unit cell along the [111] direction. As a result, there are 
three crystallographically inequivalent V sites as shown in Fig. 1. Bond lengths between V 
sites are estimated as 3.04, 2.81, 2.61 and 3.14 A for V1-V3, V2-V3, short V3-V3 (within 
trimers) and long V3-V3 (between trimers) bonds, respectively. 10 

These substantial differences among the bond lengths imply a decomposition of the whole 
pyrochlore lattice into two components; one is 'heptamer', i.e., a cluster with seven V sites 
(one V2 site and six V3 sites) connected by bold and gray lines in Fig. 1, and the other is a 
remaining VI site. This decomposition and the heptamer formation may provide an alternative 
scenario to the previous charge-ordering one, which we call the heptamer scenario hereafter. 
The formation of the heptamer itself naturally explains the doubling of the unit cell along 
the [111] direction. The peculiar temperature dependence of the susceptibility may also be 
compromised if one assumes spin-singlet states in the heptamers and local moments at each 
VI sites; the former leads to a sudden drop at T c and the latter gives rise to the Curie-like 
component below T c . In fact, it was pointed out that the latter Curie-like component is well 
fitted by assuming an S = 1 local moment at each VI site. 10 

From a theoretical point of view, two important issues must be addressed to understand 
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the physical properties of AIV2O4 on the basis of the heptamer scenario, (i) Why are the 
heptamers stabilized in this frustrated pyrochlore system? (ii) How does the spin-singlet state 
emerge in each heptamer? To answer the first question, it is necessary to consider charge, 
spin and orbital degrees of freedom in the entire pyrochlore lattice, which is very complicated. 
The second question is also nontrivial. There are in total 18 d electrons in each heptamer 
when we assign 2 d electrons to each VI site to form the S = 1 moment for the Curie-like 
contribution. It is not trivial whether we obtain a spin-singlet state with 18 electrons on 
the seven-site cluster when we take account of the charge, spin and orbital states as well as 
electron correlations. 

In this paper, we will investigate the second problem, that is, the mechanism of spin-singlet 
formation, by studying the electronic state of one heptamer explicitly. We will discuss the 
effects of orbital-dependent transfer integrals, electronic correlations and the crystal field, and 
clarify the parameter region where the spin-singlet ground state is obtained. This systematic 
study provides a simple physical picture of the spin-singlet formation in the heptamer. We 
believe that the stabilization mechanism of the spin-singlet state may give a hint for the first 
question as well. 

This paper is organized as follows. In §2, we introduce a heptamer Hamiltonian on the 
basis of a multiorbital Hubbard model. Two different singlet states in the noninteracting 
case are discussed in §3. In §4, we show our numerical results for the electronic state of the 
heptamer Hamiltonian. Section 5 is devoted for the summary and concluding remarks. 

2. Heptamer Hamiltonian 

To consider the electronic state of one heptamer, we start from a multiorbital Hubbard 
model with three-fold t2 g orbital degeneracy. The Hamiltonian consists of three terms as 

Wf u ll = Hkin + Hint + Htrig, (1) 

where Wkin is the kinetic term determined by electron hoppings between V sites, Ti^ rep- 
resents the on-site Coulomb interactions, and Htrig is for the crystal- field splitting of tig 
levels induced by the trigonal distortion of VOg octahedra. The Hamiltonian is defined on the 
seven-site cluster in Fig. 2(a), and contains 18 electrons in total as mentioned in §1. 
The first term Wkin is given in the form 

H "» = - £ £ £ t]/ (4 7T C, VT + H.C.) , (2) 
(ij) 7,7' r 

where i,j and r(=|,J,) are site and spin indices, respectively, and 7,7' = 1 (d xy ), 2 (d yz ) 
and 3 (d zx ) are orbital indices. The site indices take the values from to 6 as shown in 
Fig. 2(a), and the summation (ij) is taken over the nearest-neighbor sites. (At this stage, 
we neglect the difference between the bond lengths of V2-V3 bonds and V3-V3 bonds in the 
heptamer.) The heptamer originally consists of the edge-sharing network of VOq octahedra 
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Fig. 2. (a) One heptamer. The site numbering used in our model is shown, (b) Edge-sharing network 
of VC>6 octahedra for the heptamer. Oxygen atoms are at the corners of octahedra. Transfer 
integrals between tig orbitals of (c) a-, (d) ~k\- and (e) 7r2-type transfer integrals, respectively. 

as shown in Fig. 2(b). For the edge-sharing configuration, the d-d transfer integrals tJJ have 
three different elements as shown in Figs. 2(c)-(e). Here, we call them the a-, 7Q- and 7r2-type 
transfer integrals. 

The Coulomb interaction term TCmt consists of the four components as 

Hint = 7~iu + 7~tu' + + Hj', (3) 

where 

H u = u YY n if1 n nl ' ( 4 ) 

i 7 

HU> = U'y^y~] THjUyyl, (5) 

i 7>7' 

Hj = 3 Y Y Y 4yT C i7'T4yV c Hr' , ( 6 ) 
i 7>7' t,t' 

i 77^7' 
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Here, Tijj, TCjj', 7~Lj and T~tj' denote the intra- and inter-orbital Coulomb interactions, the ex- 
change interaction and the pair-hopping term, respectively. The density operators are defined 
by n ilT = c\ 1T Ci JT and n i7 = J2 T n iir- 

The third term in eq. (1), TCtrig, represents the effect of the trigonal distortion of YOq 
octahedra. The lattice structure determined by experiments indicates that the VC>6 octahe- 
dron including the V2 site is substantially compressed in the [111] direction. 10 The trigonal 
distortion leads to the crystal-field splitting of three-fold tig levels at the V2 site into the 
low-energy a\ g singlet and the high-energy e g doublet. This crystal-field effect is described by 



The compression of VOg corresponds to a positive D, and the energy splitting between the a\ g 
singlet and the e g doublet is 3D/2 in this definition. We neglect small distortions of octahedra 
including V3 sites for simplicity. 

The multiorbital Hubbard model given by eq. (1) is too complicated to fully handle its 
electronic state, but we can make a simplification on the basis of the lattice structure de- 
termined by the experiment. The key observation is that the length of the V3-V3 bonds is 
considerably shorter than that of V2-V3 bonds. In fact, the V3-V3 bonds in the heptamer are 
the shortest ones in the original pyrochlore system, constituting the trimers in the Kagome 
layers. Hence, we expect that the cr-type transfer integrals for these V3-V3 bonds are the 
most dominant among various contributions in tJJ . From this observation, we assume that 
the cr-type bonding state for each V3-V3 bond has a sufficiently low energy and is fully oc- 
cupied by two d electrons per bonding state, forming a singlet. For example, the site i = 1 is 
on two V3-V3 bonds = (1,2) and (1,3) which lie in the yz and zx planes, respectively, 
and therefore, the d yz and d zx orbitals constitute the bonding levels. This is indeed the case 
in the noninteracting limit 7i\ n t = 0, and is supposed to be valid in the weak-correlation 
regime. The interorbital interactions at V3 sites will be taken into account within the Hartree 
approximation [see eq. (14)]. 

The assumption of the dominant cr-type V3-V3 bonds allows us to simplify the model 
defined by eq. (1) in the following points, (i) Since there are six V3-V3 bonds in the heptamer, 
6 x 2 = 12 electrons occupy the cr-type bonding states. These 12 electrons become inactive 
on the cr-type V3-V3 bonds. Because we have in a total of 18 electrons in the model given 
by eq. (1), we will consider the remaining 6 electrons in the following calculations. Hereafter, 
we call the orbitals constituting the cr-type V3-V3 bonds inactive orbitals and the remaining 
orbitals active orbitals. (ii) All 7T2-type transfer integrals [Fig. 2(e)] can be neglected, because 
orbitals which are connected to active orbitals through 7r2-type bonds are always inactive, 
(iii) 7i"i-type transfer integrals for V2-V3 bonds can also be neglected on the basis of a reason 
similar to that in item (ii). As a consequence of these simplifications, among the transfer 




(8) 



7>7' r 
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integrals, only two contributions are relevant; one is the cr-type transfer integral for V2-V3 
bonds, which is denoted by iy2-V3' an< ^ the other is the 7Ti-type one for V3-V3 bonds, iv3-V3- 
Consequently, on the basis of the above assumption, the Hamiltonian given by eq. (1) is 
simplified into the form 

T~L = ^kin + Wint + Wtrig- (9) 

The kinetic term for the remaining 6 electrons, Hkin, is given by 

rC, -1/V2-V3 i -T/V3-V3 / ln \ 

where T~C^ V3 denotes the cr-type hoppings between the V2 site (i = 0) and V3 sites (i = 1 — 6), 
and H^fa V3 denotes the 7Ti-type hoppings between the nearest-neighbor V3 sites: 

6 

T^kin V3 = _i V2-V3 E E {^OaiT^aiT + H.C^ , (11) 

1=1 T 

= -*V3-V3 E E(4^ + H.c). (12) 

fo')GV3 r 

Here, a» denotes the active orbital at the i-th V3 sites [for instance, a± = 1 (d X j/)]- Note 
that in eq. (11), the orbital a» constitutes the cr-type bond between 0-th and i-th sites, while 
the orbitals aj and aj at neighboring V3 sites constitute the 7Ti-type bond in eq. (12) where 
the summation is taken over the nearest-neighbor V3 sites. The Coulomb interaction term in 
eq. (9) is written as 

/tint — H-int + n inf 

The former term is for the V2 site, given by the terms in eq. (3) for i = only. The latter is 
for the V3 sites, which is reduced into 

6 6 

Wtat = U E nia^n iad + (2U' - J) E "i* > ( 14 ) 

i=i i=i 

where «j denotes the active orbital as in eqs. (11) and (12). Here, the second term is the 
Hartree potential from electrons in the cr-type bonding states assumed for the V3-V3 bonds. 

3. Noninteracting Case: Two Different Scenarios for Singlet Formation 

Before going into the numerical study of the electronic state of the heptamer model given 
by eq. (9), here we consider the noninteracting case, i.e., Wj nt = (U = U' = J = J' = 0). In 
this case, we have three parameters ty2-V3' ^V3-V3 ano - 

There are two important cases in this parameter space to consider the singlet formation 
in the heptamer. One is ty2-V3 ^ ano - ^V3-V3 = D = 0. The cr-type transfer integrals in an 
orbital sector (say d xy ) have finite matrix elements on one straight V3-V2-V3 bond lying in 
the corresponding plane (in the xy plane). Hence, in this case, three-site clusters are formed on 
straight V3-V2-V3 bonds in each orbital sector, for instance, a cluster with the sites i = 1, 0, 4 
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Fig. 3. Schematic pictures of two different singlet states in the noninteracting case, (a) cr-singlet state 
[eq. (15)] and (b) 7r-singlet state [eq. (16)]. Blue, red and green circles represent d xy , d yz and d zx 
orbitals at each V site, respectively. Singlet bonds are shown by bold lines. In (b), at the V2 site, 
two electrons occupy the a± g singlet state due to the trigonal splitting. 



in the d xy orbital sector. In each cluster, the ground state is given by the bonding state of three 
sites, being singlet as shown in Fig. 3(a). The ground-state wave function for the heptamer is 
given by 



3 1 

nn 2 (4*^+^4^+ 4 



(i+3)aiT 



10), 



(15) 



t i=l 



where |0) is the vacuum state of the model in eq. (9). We call this the cr-singlet state. Note 
that the state in eq. (15) is obtained for both positive and negative tv2-V3- 

The other important case is iy2-V3 = ^ anc ^ ^V3-V3'^ ^ 0. In this case, the heptamer is 
separated into the isolated V2 site and two trimers with V3 sites. In each trimer, the ground 
state is given by the bonding state due to the 7ri-type transfer integrals, iy3-V3- At the V2 
site, for a positive D, the ground state is given by a doubly occupied a\ g singlet as shown in 
Fig. 3(b). Hence, the wave function for the heptamer is expressed as 

3 J 



3 J 3 t 6 t 

1 ; il h vsti vs vs 



|0). 



(16) 



t 7=1 ' ~ 8=1 ' ~ 1=4 

This state is schematically shown in Fig. 3(b). We call this the 7r-singlet state. Note that if 
*V3-V3 or D is negative, the ground state no longer becomes singlet and has some degeneracy. 

When we switch on the interactions, it is nontrivial that the ground state is singlet. 
Nevertheless, as we will see in the numerical results in the next section, there appear two 
different singlet regions; one includes the cr-singlet limit and the other includes the 7r-singlet 



8/17 



J. Phys. Soc. Jpn. Full Paper 

one. In other words, in each region, the wave function is adiabatically connected with either 
|cr) or \tt). Besides these two singlet regions, there appears also a doublet nonmagnetic ground 
state. Thus, our questions are twofold, (i) For realistic values of parameters, does the ground 
state of the heptamer Hamiltonian become singlet? (ii) If that is the case, which category 
does the singlet state belong to, the \a)- or |-7r)-like state? 

4. Effects of Electron Correlation: Exact Diagonalization Study 

4-1 Method and parameters 

We investigate the electronic state of the model in eq. (9) including the interaction terms 
by exact diagonalization (ED). We use the Householder method for the ED calculations. 

For the d-d transfer integrals, the Slater-Koster scheme 11 gives the estimates ty 2 V3 ~ 
0.5 eV and £y3-V3 — —0.26 eV on the basis of using the bond lengths in experiments. Hereafter, 
we set ty2-V3 = 1 as an energy unit, and change *V3-V3 f rom —2 to 2 and D from to 2 to 
study overall features in the parameter space, which is supposed to include realistic values for 
the present tig system. 

Concerning the interaction parameters, no experimental estimate for AIV2O4 is available 
yet. An estimate was given for a related spinel vanadate L1V2O4; 12 the optical measurement 
suggests U ^ 2eV. In the following, we present the results for U = — 8 to show the tendency 
with changing U. We set J = 0.1U, and retain the relations U = U' + 2 J and J = J'. 13 
4-2 Results 

First, we study the magnetism and the degeneracy of the ground state of the heptamer 
Hamiltonian given by eq. (9). Figure 4 shows the results in the iy 3 V3 -.D plane for different 
values of U. Nonmagnetic ground states are obtained in three regions in the present parameter 
space; two singlet regions without any degeneracy (red and blue) and a doublet region (green). 
Two singlet regions are well separated, and each region includes one of the limits in the 
noninteracting case discussed in §3; the red region extends from the cr-singlet limit (t V 3-V3 = 
D = U = 0), the blue region continues to the 7r-singlet limit (iy2-V3 = U = and ty3-V3' D > 
0). The doublet region emerges only for finite values of U. Both singlet regions shrink as 
U increases, but remain well for intermediate or rather large values of U, while the doublet 
region spreads as U increases. 

We focus on two singlet regions, one of which is believe to be relevant to understanding the 
experiments in AIV2O4 as discussed below. We will discuss the doublet region later. Although 
these two singlet regions are separated in the parameter space, the ground states for two 
regions show the same symmetry. The symmetry of the ground state of the model given by 
eq. (9) belongs to one of the irreducible representations Ai, A2 and E because the Hamiltonian 
belongs to the D3 point group. 15 We calculate (4>o\Cs\4>o) an d (^olC^o), where C3 and C 2 
are operators for the symmetric rotation of a threefold symmetric axis by 3ir/2 and that of a 
twofold symmetric axis by n, respectively. Here, \4>o) is the ground-state wave vector obtained 
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by ED calculations. We found that (</>o|C3|0o) = {4>o\C' 2 \(j)o) = 1 for all the ED solutions 
without degeneracy. Namely, all the ground states in both separated singlet regions belong to 
the A\ symmetry group. 

To characterize two singlet regions, we calculate overlaps of the ground-state wave function 
with \a) and \ir) in eqs. (15) and (16), respectively. We plot the square of overlaps in Fig. 5 
for typical sets of parameters. In the red region in Fig. 4, the square of the overlap with \a), 
(<r|(/>o) 2 , is substantial whereas (vr|0o) 2 is negligibly small (< 0.01). Since (<t|0o) 2 3> (ir\4>o) 2 , 
we call this red region the a-singlet region. On the other hand, in the blue region in Fig. 4, 
(tt\4>o) 2 has a substantial value compared with (a\(j)o) 2 < 0.01, and therefore we call the blue 
region the 7r-singlet region. Hence, these two singlet regions are well distinguished by the values 
of overlaps, and appear to be separated within the parameter space that we investigated. 

We also calculate the charge disproportionation between the V2 site and the V3 sites in 
the ground state. Figure 6 shows the electron density at the V2 site; 



Note that the electron density at the V3 site is given by (18 — A^2)/6 for singlet states with 
A\ symmetry. In the a-singlet region (red region in Fig. 4), as U increases, N 2 continuously 
decreases from the noninteracting value 3. On the contrary, in the 7r-singlet region (blue region 
in Fig. 4), N 2 is almost unchanged and remains at the noninteracting value 2. This difference 
comes from the nature of the singlet states. In the former, since the dominant singlet bonds 
are on the straight V3-V2-V3 clusters, electron charges are easily redistributed between V2 
and V3 sites. On the other hand, in the latter 7r-singlet region, the dominant singlet bonds 
are confined within V3 trimers, and the bonds between V2 and V3 sites are relatively weak. 
There, the singlet energy gain within the V3 trimers as well as the crystal field D at the V2 
site prevents the charge redistribution between V2 and V3 sites. 

The results in Figs. 4 and 5 show that for the estimates of transfer integrals by using the 
Slater-Koster scheme, the ground state of the model given by eq. (9) is singlet and belongs 
to the cr-singlet region where the ground-state wave function is close to \a) in eq. (15). We 
note that the 7r-singlet region is only limited for large values of iv2-V3 an( ^ ^ compared with 
ty 2 _Y 3 , and this is not realistic in the present t 2g electron system. Therefore, we conclude that 
the singlet ground state is realized in the heptamer, and that the basic mechanism of singlet 
formation is understood on the basis of the cr-type bonding of t 2g orbitals. 

We calculate the spin gap A for a comparison with experimental data, which we define 
here as the energy gap from the ground state to the lowest magnetic excited state. Note that 
the spin gap becomes zero when the ground state is magnetic by definition. The result is shown 
in Fig. 7 for several Slater-Koster parameters. Although a quantitative comparison with the 
experimental estimate is beyond the scope of the present study since interactions among the 




(17) 
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heptamers are completely neglected in the present theory, 16 this finite spin gap gives rise to the 
suppression of magnetic susceptibility below T c observed in experiments. Therefore, our results 
on the gapped, spin-singlet ground state in the heptamer model support the heptamer picture, 
in which the peculiar temperature dependence of the magnetic susceptibility is explained as 
the summation of the heptamer-singlet contribution and the S = 1 Curie-like contribution. 10 
Finally, let us make a remark on the doublet state (green region) in Fig. 4. The doublet 
region develops as U increases, and includes the parameter regime given by the Slater-Koster 
scheme for large values of U > 8. From the symmetry analysis, we find that the doublet state 
belongs to E symmetry. Since this state is also nonmagnetic, it might be possible to explain the 
peculiar temperature dependence of the magnetic susceptibility. Nevertheless, we believe that 
the u-type singlet ground state is relevant to AIV2O4 for the following reasons. One is that 
the degeneracy remaining in the doublet state should be lifted by the generalized Jahn- Teller 
mechanism or the heptamer-heptamer interaction when we allow a further lattice distortion 
of the heptamer; however, such lowering of the symmetry is not observed in experiments. 
The other reason is that our present analyses are justified in the weak correlation regime as 
discussed in §2, and hence some effects beyond our assumptions may alter the results in the 
relatively large U region. 

5. Summary and Concluding Remarks 

We have investigated the origin of spin-singlet formation in seven-site clusters, heptamers, 
which emerge in the spinel oxide AIV2O4 as a result of the structural change with the trimer 
formation in Kagome layers of the V pyrochlore lattice. We have derived an effective model to 
describe the charge, spin and orbital degrees of freedom in one heptamer from the Hubbard 
model with t2 g orbital degeneracy, and investigated its electronic and magnetic states using 
numerical exact diagonalization. Our results indicate that in the realistic parameter region 
for AIV2O4, the ground state of the heptamer becomes spin-singlet, and that not the trimer 
formation in the Kagome layers but the heptamer formation plays an essential role in the 
spin-singlet formation in this system. We conclude that the singlet state is stabilized by 
dominant d-d transfer integrals of the cr-type which form bonding states within the heptamer. 
The heptamer formation can give a comprehensive understanding of the magnetic properties 
observed in experiments; a suppression due to the singlet gap formation and the Curie-like 
behavior at lower temperatures coming from 5=1 moments outside the heptamers. 

How can we distinguish experimentally the three different origins of spin-singlet formation, 
i.e., two different singlet states (a- and 7r-type) and one doublet state? We can distinguish 
between singlet and doublet states by the entropy at high temperatures estimated from the 
electronic specific heat. Then, it is possible to distinguish between a- and 7r-type singlet states 
which have the same symmetry. One possible way is to detect the directional dependence of 
the spin excitation. The dominant singlet bonds lie on the intertrimer V3-V2-V3 bonds in the 
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a case and within the V3 trimers on Kagome planes in the n case. This difference may be 
detected using the Raman scattering technique. Another possibility is to measure the local 
density of d electrons at the V2 site. Our results in Fig. 6(f) indicate that the local density 
becomes ~ 2.5 — 3.0 in the cr-singlet region, while it is ~ 2.0 in the 7r-singlet one. 

Let us comment on the properties of S = 1 moments outside the heptamers. In our hep- 
tamer scenario, below T c , the magnetic moments that remain active are the S = 1 spins at VI 
sites only, and the spins at V2 and V3 sites are magnetically inactive, because of the singlet 
formation in heptamers. Since a VI site is connected to the V2 sites (see Fig. 1), the VI spins 
are almost magnetically isolated. At a much lower temperature than T c , another anomaly 
appears at T ~ 5K in the magnetic susceptibility, which is ascribed to spin-glass transi- 
tion. 8 This indicates that in this low-temperature region, S = 1 spins begin to interact with 
each other. Nearest-neighbor VI pairs correspond to the third-nearest-neighbor sites across 
hexagons in the original pyrochlore lattice, forming a two-dimensional triangular lattice in the 
[111] planes. Hence, the S = 1 moments at VI sites constitute a weakly coupled triangular 
spin system at low temperatures. The sign of exchange interactions is not trivial because of 
complicated exchange paths, but in any case, it is believed that Heisenberg spin systems on a 
weakly coupled triangular lattice show a magnetic order, not a spin-glass behavior. Therefore, 
our heptamer scenario suggests that the spin-glass behavior comes from an extrinsic effect, 
such as randomness. 

Our present study was conducted on the basis of the lattice structure below the structural 
transition temperature, and therefore, at this stage, we cannot answer why and how the lattice 
distortion including the formation of heptamers is stabilized in the frustrated pyrochlore lattice 
structure. There are many similar examples of such clusters in frustrated systems, such as 
hexamers in ZnCr204, 17 octamers in CuIr2S4, 5 ' 18 trimers in L1VO2, 19 and dodecamers in a 
double-exchange spin-ice model on a Kagome lattice. 20 These facts suggest the ubiquitous 
role of cluster formations in lifting the degeneracy inherent in frustrated systems, although 
the driving force may depend on the details of systems. In the present system AIV2O4, our 
results imply that the stability of the heptamer singlet state itself may play some role in 
the mechanism of self-organizing the clusters. However, to identify a dominant player, it is 
necessary to handle charge, spin, orbital and lattice degrees of freedom and to compare the 
instabilities to many possible ordered states on the frustrated pyrochlore lattice. This is left 
for future study. 
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